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Abstract 


A  finite-element  analysis  is  carried  out  for  small -amplitude  free 
vibration  of  laminated,  anisotropic,  rectangular  plates  having  arbitrary 
boundary  conditions,  finite  thickness-shear  moduli,  rotatory  Inertia,  and 
bimodulus  action  (different  elastic  properties  depending  upon  whether  the 
fiber-direction  strain  is  tensile  or  compressive).  The  element  has  five 
degrees  of  freedom,  three  displacements  and  two  slope  functions,  per  node. 

An  exact  closed-form  solution  Is  also  presented  for  the  special  case  of 
freely  supported  single-layer  orthotropic  and  two-layer,  cross-ply  plates. 
This  provides  benchmarks  to  evaluate  the  validity  of  the  finite-element  anal¬ 
ysis.  Both  solutions  are  compared  with  numerical  results  existing  in  the 
literature  for  special  cases  (all  for  ordinary,  not  bimodulus,  materials)  and 
good  agreement  is  obtained. 


Presently,  Professor  of  Engineering  Science  and  Mechanics,  Virginia 
Polytechnic  Institute  and  State  University,  Blacksburg,  Virginia. 
^Presently,  Structures  Engineer,  Lear  Fan  Corp.,  Reno,  Nevada. 
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Introduction 

Structural  uses  have  been  increasing  for  laminates  consisting  of  multi¬ 
ple  layers  of  fiber-reinforced  composite  materials.  Consequently,  there  is 
an  increasing  need  for  more  realistic  mathematical  modeling  of  the  material 
behavior  for  incorporation  into  static  and  dynamic  structural  analyses. 

Certain  fiber-reinforced  materials  have  been  found  experimentally  to  exhibit 
quite  different  elastic  behavior  depending  upon  whether  the  fiber-direction 
strain  (e^)  is  tensile  or  compressive  [1-3].  Examples  of  such  materials  are 
tire  cord-rubber,  reinforced  solid  propellants,  and  some  biological  tissues. 
Although  the  stress-strain  behavior  of  such  materials  is  actually  curvilinear, 
it  is  often  approximated  as  being  bilinear,  with  different  slopes  (elastic 
properties)  depending  upon  the  sign  of  ef.  Thus,  they  are  called  bimodulus 
composite  materials. 

The  limited  number  of  previous  analyses  of  bimodulus-material  plates 
were  reviewed  in  [4-6],  and  all  were  limited  to  static  analyses.  Thus,  it  is 
believed  that  the  present  work  is  the  first  vibrational  analysis  of  such 
plates.  The  present  work  Is  not  limited  to  just  thin  plates  of  isotropic  bi- 
modulus  material,  rather  it  is  applicable  to  moderately  thick  plates  laminated 
of  orthotropic  bimodulus  material.  Two  formulations  are  presented  and  solved: 
one  is  a  mixed  finite-element  formulation  with  five  degrees  of  freedom  per 
node,  and  the  other  is  an  exact  closed-form  solution. 

Governing  Equations 

Mindlin's  linear  dynamic  theory  [7]  of  moderately  thick  plates  was  first 
extended  to  plates  laminated  of  ordinary  (not  bimodulus)  monoclinic  elastic 
material  by  Yang,  Norris,  and  Stavsky  (YNS)[8].  Later,  Wang  and  Chou  [9] 
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showed  that  a  slightly  different  version  of  the  YNS  theory,  presented  by 
Whitney  and  Pagano  [10],  Is  more  accurate  than  the  original  version  [8]. 

Here,  this  class  of  theory  is  extended  to  bimodulus-material  laminates. 

The  origin  of  a  Cartesian  coordinate  system  is  taken  to  be  In  the  mid¬ 
plane  (xy  plane)  of  the  plate  with  the  z  axis  being  normal  to  this  plane  and 
directed  positive  downward. 

Using  the  fiber- governed  symmetric  material  model  Introduced  In  [11], 
we  take  the  generalized  Hooke's  law  for  the  in-plane  action  In  each  layer  (t) 
to  be  of  the  following  bimodular  form: 


Here  the  stresses  (o  ,c  ,t  )  and  engineering  strains  (ev,eu,yv  )  are  denoted 

a  y  a  y  x  y  xy 

In  the  usual  fashion,  and  the  Q's  are  the  plane-stress-reduced  stiffnesses 
(symmetric  array).  The  first  two  subscripts  of  the  Q's  are  those  classically 
used  in  anisotropic  elasticity  [12]  and  composite-material  mechanics  [13]. 
Here  the  third  subscript  (k)  refers  to  the  sign  of  the  fiber-direction  strain 
( k»l  for  tension  and  k=2  for  compression),  and  i  refers  to  the  layer  number 
Ual,2,...,n,  where  n  is  the  total  number  of  layers). 

It  is  assumed  that  the  thickness-shear  behavior  is  unaffected  by  bi¬ 
modular  action,  thus 


The  stress  and  moment  resultants,  each  per  unit  length,  are  expressed 


in  terms  of  stresses  as 


rh/2 

(Nx,Ny.Nxy*Qx»Qy) 

i-h/2  (ox,°y,Txy*Tx2*Ty2)  dZ 

(3) 

rh/2 

<WTxv)z  dz 

j-h/2  x  y  xy 

(4) 

The  displacement  components,  u,  v,  and  w  in  the  x,  y,  and  z  directions, 
respectively,  can  be  expressed  in  terms  of  mid-plane  displacements  u°,  v°, 
w°  and  slope  functions  <|»x  and  ip  as: 


=  u°(x,y,t)  +  z^x(x,y,t)  ;  v  *  v°(x,y,t)  +  z*y(x,y,t) 

w  «  w(x,y,t) 


(5) 


where  t  is  time. 

The  constitutive  equations  for  an  unsymmetric  cross-ply  laminate  are: 


xy 


>  * 


xy 


Au 

a12 

0 

Bu 

b12 

0 

a12 

^22 

0 

B12 

B22 

0 

0 

0 

^66 

0 

0 

b66 

.-H 

GO 

Bl2 

0 

°12 

0 

Bl2 

^22 

0 

O12 

O22 

0 

0 

0 

Bee 

0 

0 

o66 

v0 
v  *\i 


+  ||U 

v,x  u’y 


rx,x 


ry,y 


J  v.  vx,y J 


(6) 


and 


S44  0 

0  S55 


w  +  \l> 

»y  y 

w’x  +  *x 


(7) 
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Here  differentiation  is  denoted  by  a  comma,  i.e. ,  (  )>x  =  3(  )/ax, 

and  the  extenslonal,  flexural -extenslonal  coupling,  and  flexural  stiffnesses 
of  the  laminate  are  defined  by 

h/2 

<A1J-B1J*V  ■  \_m  (8) 

1»j*l»2,6 

Also,  the  thickness-shear  stiffnesses  of  the  laminate  are  defined  by 

sti  ■  kich  <-4-5  <«> 

where  the  K*  are  the  thickness-shear  correction  coefficients,  which  can  be 
determined  by  various  approaches,  cf.  [12]. 

In  addition  to  performing  the  integrations  in  a  piecewise  manner  from 
layer  to  layer,  one  also  has  to  take  into  consideration  the  possibility  of 
different  elastic  properties  (tension  or  compression)  within  a  layer.  This 
is  explained  in  detail  for  a  two-layer  cross-ply  laminate  in  Appendix  A. 

Taking  into  account  the  coupling  and  rotatory  inertias,  one  can  write 
the  equations  of  motion  as  follows: 


N..  ..  +  N 


Pu?„  +  R*v 


x,x  xy,y  '“’tt  'rx,tt 
Nxy,x  +  "y,y  "  Pv?tt  +  *y.tt 
Qx,x  +  Qy,y  " 

Mx,x  +  Mxy.y  "  Qx  *  Ru’tt  +  r*x,tt 
Mxy,x  +  My,y  "  ^y  *  Rv’tt  +  ^y.tt 


(10) 


Here  P,  R,  and  I  are  the  normal,  coupling,  and  rotatory  inertia  co¬ 


efficients  per  unit  mid-plane  area  and  are  defined  by 
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fh/2 

(P.R.I)  -  eO.z.z2)  dz  (11) 

J-h/2 

where  p  Is  the  material  density. 

Substituting  equations  (6)  and  (7)  into  equations  (10),  we  obtain  the 
equations  of  motion.  In  operator  form,  we  have 


where  [L^]  is  a  symmetric  linear  differential  operator  matrix  with  the 
following  elements: 

Ln  5  And2  +  A66d2  -  Pd£  ;  l12  =  (A12  +  A66)dxdy  ;  L13  s  0 

Ll4  =  C(B12  +  B66)/h]dxdy  ;  L15  =  (Bu/h)d2  +  (B66/h)d2-(R/h)d2 

L22  =  A66d2  +  A22dJ-PdJ  ;  l23  =  0  ;  l24  =  (B66/h)d2  +  (B22/h)dJ  -  (R/h)d| 


L25  =  U**  »  L33  5  “  S5sdx  -  S44dy  +  Pd2 

L34  2  -  (S44/h)dy  ;  L35  =  -  (S55/h)dx 
L44  ^  (D66/h2)d2  *  (D22/h2)d2  -  (S44/h2)  -  ( I/h2)d2 


(13) 


Us  *  C(0u  +  06s)/h23dxdy 

Us  2  COn/h2)d2+(D66/h2)dy-(S5S/h2)-  (I/h2)d2  ;  dx  =  a/ax,  etc. 
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Application  to  Plate  Freely  Supported  on  all  Edges 

The  boundary  conditions  on  all  edges  are  freely  supported  (simply  sup¬ 
ported  without  in-plane  normal  restraint). 

Along  the  edges  at  x  *  0  and  x  =  a. 


w  ® 


M. 


Nx.O 


Along  the  edges  at  y  *  0  and  y  =  b. 


(14) 


w  *  *x  =  My  “  0 
u°  -  Ny  =  0 


Closed-Form  Solution 

The  governing  equations  (12)  and  the  boundary  conditions  (14)  are 
exactly  satisfied  in  closed  form  by  the  following  set  of  functions: 

u°  *  U  cos  ax  sin  By  eiu)t 

v°  *  V  sin  ax  cos  gy  ei<4it 

w  s  W  sin  ax  sin  gy  e1uit  (15) 

hij>y  *  Y  sin  ax  cos  gy  e1u)t 

*  X  cos  ax  sin  gy  e1u,t 

Here,  w  is  the  natural  frequency  associated  with  the  mode  having  axial  and 
transverse  wave  numbers  m  and  n,  and 

a  =  nnr/a  ,  B  =  mr/b  (16) 

where  a  and  b  are  plate  dimensions  In  the  x  and  y  directions,  respectively. 
Substituting  solutions  (15)  into  the  governing  equations  (12),  we 


obtain  the  following: 
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V 

'o' 

V 

0 

<  w 

*  4 

0  - 

Y 

0 

(17) 


k»£sl ,2,3, 4, 5 


where  Ckjl  is  a  5x5  symmetric  determinant  containing  the  following  elements: 

C 1 1  =  -  Ana2  -  Aggg2  +  Pu2  ;  C12  =  -  (Ai2  +  A6g)a8  ;  C 1 3  =  0 

C14  =  -  [(Bi2  +  B66)/h]a8  ;  C15  =  -  (Bn/h)a2- (B66/h)B2+(R/h)<,2 

C22  -  "  Agga2  -  A22B2  +  Pw2  »  C23  -  0 
C24  -  ~  (B66/h)a2  -  (B22/h)82  +  (R/h)w2  i  C2S  =  CiJ# 

(18) 

C33  =  -  (S55a2  +  S4482  -  Pu2)  ;  C34  =  -  (S44/h)g 

C35  =  -  (S55/h)a  ;  C44  =  -  (Dgg/h2)a2.(D22/h2)B2-(S44/h2)  +  (I/h2)a,2 

C45  =  -  [{D12 + Dgg)/h2]aB 

C55  s  -  (Dn/h2)a2  -  (Dgg/h2)s2  -  ( S 55/h2 )  +  (I/h2)w2 

The  frequency  u  can  be  determined  by  setting  lckJ  *  °* 

To  determine  the  z-posltion  of  the  fiber-direction  neutral  surface,  one 

sets 

s  e°  +  Zk^  *  0 

znf  *  "  zVKf  ' 


or 


(19) 
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Thus,  z  =  -  hU/X  and  z  =  -  hV/Y.  An  iterative  procedure  is  used  to  obtain 

IIA  iljF 

the  final  displacement  ratios  and  corresponding  frequency. 


Finite-Element  Formulation 

An  exact  closed-form  solution  to  equations  (12)  can  be  obtained  only 
under  special  conditions  of  geometry,  edge  conditions,  loadings,  and  lamination. 
Here  we  present  a  simple  finite-element  formulation  which  does  not  have  any 
limitations  (except  for  those  implied  in  the  formulation  of  the  governing 
equations) [15]. 

Suppose  that  the  region  His  subdivided  into  a  finite  number  N  of  sub- 
regions:  finite  elements,  H  (e«l,2,...,N).  Over  each  element  the  general¬ 
ized  displacements  (u°,v°,w,<|/  ,i|>  )  are  interpolated  according  to 

x  y 


o*i  o  i  *  2 

U  a  l  Uj+j  ,  V  *  r  V^^  ,  W  *  Z  w^ 

P  3  P  3 

*x  =  {  Vi  •  *y  =  *  Vi 


(20) 


where  (a=l,2,3)  is  the  interpolation  function  corresponding  to  the  i-th 
node  in  the  element.  Note  that  the  in-plane  displacements,  the  transverse 
displacement,  and  the  slope  functions  are  approximated  by  different  sets  of 
interpolation  functions.  While  this  generality  is  included  in  the  formu¬ 
lation  (to  indicate  the  fact  that  such  independent  approximations  are  possible), 

we  dispense  with  it  in  the  Interest  of  simplicity  when  the  element  is  actually 

12  3 

programmed  and  take  ^  (r=s=p).  Here  r,s,  and  p  denote  the  number 

of  degrees  of  freedom  for  each  variable.  That  is,  the  total  number  of  degrees 
of  freedom  per  element  is  2r  +  s  +  2p. 


r 
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where 


,Sn 

ij 


,?n  = 

l1J 


t ?n  = 
■  j 


,!  i\ 
i»5  J»n 

dxdy 

(1.J-1.2,. 

...r) 

1  * 3 
i,€*j,n 

dxdy 

(i=l  ,2, . . . 

,r  j  j-1 ,2, . . 

...t) 

1  2 

I  .  d)  . 

l.rj.n 

dxdy 

( i*l »2» • • . 

,r  ;  j=l ,2, . . 

i  • ,  s ) 

2  2 
i.C*j»n 

dxdy 

(i »j=l .2,. 

•  •  ,s) 

2  A3 

dxdy 

(1*1,2,... 

, s  i  j- 1 ,2,. . 

...t) 

3 

i,S9J,n 

dxdy 

(i,j=l,2,. 

. . ,  s ) 

(24) 


(C,n*0,x,y) 

XX  X  12  3 

and  =  G^.,  etc.  In  the  special  case  in  which  ^  =  $.  -  <j>^,  all  of  the 
matrices  in  equations  (24)  coincide. 

In  the  present  study,  elements  of  the  serendipity  family  are  employed 
with  the  same  interpolation  for  all  of  the  variables.  The  resulting  stiffness 
matrices  are  20  by  20  for  this  four-node  element  and  40  by  40  for  the  eight- 
node  element.  Reduced  integration  [16]  must  be  used  to  evaluate  the  matrix 
coefficients  in  equations  (23).  That  is,  if  the  four-node  rectangular  element 
is  used,  the  1x1  Gauss  rule  must  be  used  in  place  of  the  standard  2x2  Gauss 
rule  to  numerically  evaluate  the  coefficients  K... 

*  J 

Substituting  solution  (22)  into  equations  (19),  we  get 


4 ’-“V**,*  •  4‘-vV*y.y 


(25) 
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Numerical  Results 

Computations  using  the  closed-form  and  finite-element  solutions  were 
carried  out  on  an  IBM  370  computer.  Since  there  is  no  previous  analysis 
for  vibration  of  bimodulus  plates,  the  present  results  could  be  compared 
only  with  those  for  rectangular  plates  laminated  of  ordinary  materials. 
Comparisons  with  the  fundamental-frequency  results  of  Jones  [17]  for  thin 
plates  and  Fortier  and  Rossettos  [18]  for  thick  and  thin  plates  are  pre¬ 
sented  in  Tables  1  and  2.  It  can  be  seen  that  the  agreement  is  good. 

As  examples  of  some  actual  bimodulus  materials,  two  composites  used  in 
automobile  tires,  aramid-cord/rubber  and  polyester-cord/rubber,  are  selected. 
The  material  properties  used  are  listed  in  Table  3.  These  are  based  on  the 
experiments  of  Patel  et  al.  [2]  and  are  the  same  data  used  in  [6]  with  the 
addition  of  the  specific-gravity  values,  which  were  estimated  on  the  basis 
of  the  volume  fractions.  The  numerical  results  for  single- layer  0°  ortho¬ 
tropic  and  two-layer  cross-ply  plates  are  presented  in  Tables  4  and  5-6, 
respectively.  As  can  be  seen  from  these  tables,  again  the  agreement  is 
good. 

There  may  be  a  question  regarding  the  effect  of  bimodulus  action  on 
plate  stiffness  in  different  portions  of  each  cycle.  For  example.  Fig.  1 
represents  a  single-layer  bimodulus-material  plate  at  the  two  extremes  of  its 
deflection.  Figure  1(a)  depicts  the  initial  half  cycle,  during  which  the  top 
surface  is  in  compression  and  the  bottom  in  tension,  thus  causing  the  neutral 
surface  for  ex  to  be  positive  (znx  >  0),  i.e.,  below  the  plate  midplane  a 
certain  distance.  Figure  1(b)  depicts  depicts  the  latter  half  cycle,  during 
which  the  top  surface  is  now  in  tension  and  the  bottom  in  compression,  thus 
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causing  znx  to  be  negative,  i.e.,  to  fall  above  the  plate  midplane.  However, 
the  absolute  value  of  znx  is  identical  to  its  value  in  the  first  half  cycle. 
Thus,  it  can  be  concluded  that  the  the  frequency  associated  with  the  second 
half  cycle  is  identical  to  that  of  the  first  half  cycle  and  either  modal 
shape.  Fig.  1(a)  or  1(b),  will  give  the  same  computational  result  for  the 
natural  frequencies. 

Now  consider  a  two-layer  laminate  with  the  bottom  layer  (layer  fc*l) 
oriented  at  0  degrees  and  the  top  layer  (t*2)  at  90  degrees;  see  Fig.  2. 
Initially,  as  shown  in  Fig.  2(a),  the  neutral  surface  for  ex  falls  below  the 
interface,  within  the  0-degree  layer,  while  the  neutral  surface  for  ey  falls 
above  the  interface,  completely  within  the  90-degree  layer.  In  the  latter 
portion  of  the  cycle.  Fig.  2(b),  the  ex  neutral  surface  falls  outside  of 
the  0-degree  layer,  and  the  ey  neutral  surface  falls  outside  of  the  90-degree 
layer.  Thus,  compressive  properties  are  used  for  the  entire  0-degree  layer, 
and  tensile  ones  for  the  90-degree  layer. 

From  the  above  considerations  for  a  two-layer  cross-ply  laminate,  it  is 
clear  that  the  plate  stiffnesses  acting  in  the  two  portions  of  a  cycle  are 
different  and  thus  the  associated  frequencies  are  also  different,  except  in 
the  case  of  a  square  plate.  Denoting  the  frequencies  associated  with  the  two 
portions  of  a  cycle  by  u1  and  to2,  it  can  be  shown  from  the  standpoint  of 
energy  conservation  that  the  average  frequency  (<*>)  over  the  entire  cycle 
must  be  given  by 

of1  =  (l/2)(w71  +  u,;1)  (26) 

Thus,  the  computational  procedure  used  for  a  cross-ply  plate  is  to  calculate 


A 
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u)i  and  u2  associated  with  modal  shapes  shown  in  Figures  2(a)  and  2(b),  respec¬ 
tively,  and  then  to  apply  equation  (26). 

Concluding  Remarks 

A  finite  element  has  been  developed  to  analyze  the  small -deflection  free 
vibration  of  laminated,  anisotropic,  rectangular  thick  plates  of  bimodulus 
material.  The  results  obtained  agree  well  with  those  of  an  exact,  closed-form 
solution  derived  for  such  a  plate  freely  supported  on  all  four  edges.  Thus, 
it  is  concluded  that  the  element  has  been  validated  and  may  be  used  for  com¬ 
putations  involving  more  complicated  boundary  conditions. 
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APPENDIX  A 


DERIVATION  OF  THE  PLATE  STIFFNESSES  FOR  TWO-LAYER 
CROSS-PLY  LAMINATE  Ot-  BIMODULUS  MATERIAL 


In  problems  involving  laminates  comprised  of  bimodulus-material  layers, 
it  is  necessary  to  evaluate  the  integral  forms  involved  in  the  definitions  of 
the  plate  stiffnesses,  equation  (8).  The  derivation  presented  here  is  for  the 
case  of  a  two-layer  cross-ply  laminate. 

Each  layer  is  assumed  to  be  of  the  same  thickness,  h/2,  and  the  same 
orthotropic  elastic  properties  with  respect  to  the  fiber  direction.  Since 
each  layer  is  oriented  at  either  0°  or  90°  to  the  x  axis,  the  laminate  is  also 
orthotropic,  i.e.,  there  are  no  stiffnesses  with  subscripts  16  and  26.  The 
bottom  layer  Is  denoted  as  layer  1,  i.e.,  a*1  in  Q^.^,  and  occupies  the  thick¬ 
ness  space  from  2*0  to  z=h/2,  where  z  is  measured  positive  downward  from  the 
midplane.  The  top  layer  is  denoted  as  layer  2,  i.e.,  £=2,  and  occupies  the 
thickness  space  from  z*-h/2  to  z=0. 

In  the  first  case  derived  here,  it  is  assumed  that  the  upper  portion  of 
the  top  layer  U=2)  is  in  compression  (k=2  in  In  the  fiber  direction 

and  that  the  lower  portion  of  the  top  layer  is  in  tension  (k*l),  while  the 
inner  portion  of  the  bottom  layer  U*l),  from  z*0  to  2*znx,  is  in  compression 
(k=2),  while  the  outer  portion  (from  znx  to  h/2)  of  layer  1  is  in  tension 
(k-1). 

Thus,  the  general  Integral  expression  for  A^,  in  equation  (8),  may  be 
taken  as  the  sum  of  the  Integrals  for  each  of  these  regions: 
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Case  1 

<0-5  >  l  >  0,  -0.5  <  l  <  0) 

x  y 

rh/2 

A1j  *  i-h/2  Qij'kJl  ^ 

z—  c0 


f  ny  fO  f*nx  fh/2 

■  J.h/2  Q1J22  dZ  *  J2  qU12  dZ  +  J„  Q1J21  dz  +  ««11  dz  <*’» 

ny  znx 

Since  the  planar  reduced  stiffnesses  are  each  respectively  constant 
in  the  appropriate  regions,  equation  (A-l)  Integrates  to  the  following  result: 


Aij  "  (Qij22  +  Q1jll)(h/2)  +  (Qij21  "  Qijll)znx 


or 


+  (Q-ij-22  "  Qijl2)zny 


Aij  *  (]/2^Qij22  +  Qijll^  +  (Qij21  '  Qijll)Zx 


(A-2) 


+  (Qij22  "  Qijl2)Zy 


(A— 3) 


Similarly 

fh/2 

*'i  '  i-h/2  "W 


rZny  rO  /Znx  rh^2 

J.h/2  ZQW2  dZ  +  Jz  ZQ1jl2  dz  +  L  zQ1j21  dz  +  j  zQ1jll  dz  (A’4) 

ny  u  znx 


-  ^1J22  +  ^1jll^h2/8)  +  (Q1 J2!  "  Q-ljVI )(znx2/2) 

+  (Q1j22  “  Q1jl2)(zny2/2) 


(A-5) 


or 
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Bij/h2  *  (1/8) (-  Qij22  +  Qijn)  +  (Qij2i  "  Qijll)(Zx2/2^ 

+  «>1j22  -  Qijl2HZy2/2) 

Also 


dz  + 


d  z  + 


dz 


=  (Q1j22  +  Q1jn)(h3/24)  +  (Qij21  ■  Qijll^Znx3/3) 

+  (Q1 j22  ‘  Qijl2)(zny3/3* 
or 

D1  j/b3  '  C!/24-) (Q1J22  +  Qijll)  +  (Q1<j21  -  Q^nK^3^) 

+  UlJ22  *  W(V/3) 


(A-6) 


(A-7) 


(A-8) 


(A-9) 


Similarly 
Case  2 

(  -0.5  <  Zx  <  0,  0.5  >  Zy  >  0) 

A1  j/h  *  (Q1jll  +  Q1j22)/2+(Qij22"<31jl2,Zx+^1j21  “Q1jll,Zy 

B1j/h2  -  *Q1jll  ■Q1j22^/8+(Q1j22'Q1jl2^Zx2/2)  +  ^1j2TQ1jll)(Zy2/2)  (A"10) 

°1j/h3  *  (91jll  +  Q1j22)/24+  <Q1 j22  *  91jl2)(Zx3/3)  +  (Q1j21  “  Q1jll )(Zy3/3) 
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Case  3 

(0.5  >  Zx  >  0,  0.5  >  Zy  >  0) 


Aij/h  *  (Q1jll  *  Qij22)/2  +  (Qij21  "  Qijll)Zx 
B1j/h2  =  (Qi jll  "  Qij22)/8  +  (Qij21  '  Qijll)(Zx2/2) 
D1j/h3  ’  (Qijl1  +  Q1j22)/24  +  (Qij21  ‘  Qijll>(Zx3/3) 


(A-11) 


Case  4 

(-0.5  <  ZY  <  0,  -0.5  <  Z  <  0) 

a  y 

Aij/h  =  (Qijl1  +  Qij22)/2  +  (Qij22  ’  Qijl2)Zy 

B1j/h2  =  «14n  -  Qlj22>/8  +  «1J22  •  Qijl2^Zy2/2>  (A‘12> 

D1j/h3  =  (Qijll  +  Qij22)/24  +  (Qij22  "  Q1jl2)(Zy3/3) 

In  the  presence  of  excessively  high  In-plane  loads,  such  as  those  due  to 
excessive  heating  of  the  midplane  or  due  to  large  deflections,  the  neutral 
surfaces  can  go  outside  of  the  thickness  of  the  laminate  and,  thus,  make  it 
act  as  it  were  homogeneous.  However,  this  does  not  occur  for  small -deflection 
free  vibrations  and  thus  the  equations  for  these  cases  are  not  presented 
here. 


Single  0°  Layer 

rnfrnrtr 

A1j/h  *  (Q1jll  +  Q1j21)/2  +  {Q1j21  "  Q1jli)Zx 

Bij/h*  -  CQijii  -  Quzi)^8  ♦  «1J21  “  ^un^x^2 

Dij/h3  *  (Q1jll  +  Q1j21)/24  +  <Q1j21  *  Q1jll)Zx3/3 


(A-13) 
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Table  1.  Comparison  of  fundamental  natural  frequencies  (m=n=l)  of  rectangular 
antisymmetric  cross-ply  plates  at  different  aspect  ratios  and  thick¬ 
nesses  (En/E22a40,  Gl2/E22=Gi3/E22=G23/E22*0.5»  v12*0>25»  K^*K|*5/6) 


Aspect 

Dimensionless  frequency 

(^b/iOV/Dzz)** 

Ratio 

a/b 

Thin-plate 

b/h- 50 

b/h=10 

theory  [17] 

C.F. 

F.E. 

C.F. 

F.E. 

0.5 

2.24 

2.400 

2.421 

1.942 

1.946 

1.0 

0.865 

0.858 

0.877 

0.794 

0.799 

1.5 

0.65 

0.656 

0.668 

0.612 

0.615 

2.0 

0.606 

0.604 

0.617 

0.565 

0.569 

2.5 

0.59 

0.590 

0.599 

0.548 

0.552 

3.0 

0.580 

0.578 

0.591 

0.541 

0.544 

★ 

C.F.  denotes  the  closed-form  solution  and 
element  solution. 

F.E.  denotes  the  finite- 

Table  2.  Comparison  of  fundamental  natural  frequencies  of  square  anti- 

symnetric  cross-ply  plates  at  different  thicknesses  (E11/E22=40, 
612/^22*^1 3/E22*l  *  G23/E22*0.5,  Vi2*0.25,  K£-K|-5/6) 


b/h 


Dimensionless  frequency  u>b2(P/E22h3),s 
Fortier  &  Rossettos  [18]  C.F.  F.E. 


10 

SO 


10.80 

11.65 


11.11 

11.82 


11.15 

12.06 
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Table  3-  Material  properties  for  two  tire-cord/rubber, 
unidirectional,  bimodulus  composite  materials3 


Aramid-Rubber 

Polyester-Rubber 

k=l  k=2 

kal 

k*2 

Longitudinal  Young's  modulus,  GPa 

Transverse  Young's  modulus,  GPa 

Major  Poisson's  ratio,  dimensionless6 
Longitudinal -transverse  shear  modulus,  GPa3 
Transverse- thickness  shear  modulus,  GPa 
Specific  gravity,  dimensionless 

3.58  0.0120 

0.00909  0.0120 

0.416  0.205 

0.00370  0.00370 

0.00290  0.00499 

0.970 

0.617  0.0369 

0.00800  0.0106 

0.475  0.185 

0.00262  0.00267 

0.00233  0.00475 

1.00 

aFiber-direction  tension  is  denoted  by  k*l,  and  fiber-direction  compression 
by  k=2. 

T. 

It  is  assumed  that  the  minor  Poisson's  ratio  is  given  by  the  reciprocal 
relation. 

3It  is  assumed  that  the  longitudinal-thickness  shear  modulus  is  equal  to 
this  one. 
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Table  4.  Dimensionless  fiber-direction  neutral -surface  locations 
and  fundamental  frequencies  for  single-layer  0°  ortho¬ 
tropic  plates  having  b/h=10  by  two  methods  (closed  form 
and  finite  element) 


Aspect 

Ratio 

a/b 


Zx  =  znx^ 


tl>b2(P/E|2h3)Js 
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Table  5.  Dimensionless  neutral -surface  locations  in  the  first  and  second  portions 
of  a  cycle  for  two-layer,  cross-ply  plates  having  b/h=10  by  closed-form 
and  finite-element  methods* 


a/b 

z  W 
y 

Z  (2> 

X 

z<2> 

y 

C.F. 

F.E. 

C.F. 

F.E. 

C.F. 

F.E. 

C.F. 

F.E. 

Aramid-Rubber: 

0.5 

0.4457 

0.4458 

-0.0648 

-0.0660 

-0.0171 

-0.0170 

0.4247 

0.4257 

0.6 

0.4446 

0.4446 

-0.0563 

-0.0554 

-0.0206 

-0.0205 

0.4303 

0.4309 

0.7 

0.4434 

0.4436 

-0.0490 

-0.0491 

-0.0240 

-0.0238 

0.4338 

0.4344 

0.8 

0.4421 

0.4421 

-0.0432 

-0.0429 

-0.0275 

-0.0274 

0.4363 

0.4365 

0.9 

0.4408 

0.4408 

-0.0385 

-0.0386 

-0.0311 

-0.0306 

0.4381 

0.4379 

1.0 

0.4394 

0.4394 

-0.0347 

-0.0344 

-0.0347 

-0.0346 

0.4394 

0.4394 

1.2 

0.4366 

0.4366 

-0.0293 

-0.0289 

-0.0424 

-0.0416 

0.4412 

0.4415 

1.4 

0.4335 

0.4337 

-0.0250 

-0.0249 

-0.0494 

-0.0497 

0.4423 

0.4426 

1.6 

0.4301 

0.4302 

-0.0218 

-0.0217 

-0.0565 

-0.0559 

0.4423 

0.4433 

1.8 

0.4265 

0.4264 

-0.0193 

-0.0193 

-0.0635 

-0.0662 

0.4437 

0.4437 

2.0 

0.4228 

0.4237 

-0.0174 

-0.0175 

-0.0705 

-0.0700 

0.4437 

0.4442 

Pol yester-Rubber : 

0.5 

0.3687 

0.3691 

-0.1335 

-0.1295 

-0.0830 

-0.0825 

0.3569 

0.357 

0.6 

0.3675 

0.3677 

-0.1213 

-0.1203 

-0.0844 

-0.0844 

0.3588 

0.359 

0.7 

0.3664 

0.3663 

-0.1119 

-0.1113 

-0.0868 

-0.0868 

0.3603 

0.360 

0.8 

0.3653 

0.3653 

-0.1050 

-0.1049 

-0.0895 

-0.0895 

0.3615 

0.362 

0.9 

0.3642 

0.3641 

-0.0999 

-0.0999 

-0.0926 

-0.0927 

0.3615 

0.362 

1.0 

0.3632 

0.3633 

-0.0960 

-0.0960 

-0.0959 

-0.0959 

0.3631 

0.363 

1.2 

0.3611 

0.3611 

-0.0906 

-0.0905 

-0.1033 

-0.103 

0.3631 

0.364 

1.4 

0.3589 

0.3596 

-0.0870  . 

-0.0870 

-0.1115 

-0.112 

0.3648 

0.365 

1.6 

0.3565 

0.3564 

-0.0846 

-0.0844 

-0.1202 

-0.120 

0.3648 

0.365 

1.8 

0.3540 

0.3538 

-0.0829 

-0.0829 

-0.1294 

-0.130 

0.3648 

0.366 

2.0 

0.3514 

0.3513 

-0.0817 

-0.0817 

-0.1389 

-0.139 

0.3660 

0.366 

*  (11 

Here  '  =  znx/h  for  the  first  portion  of  a  cycle,  etc. 
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Table  6.  Dimensionless  fundamental  frequencies  in  the  first  partial  cycle, 
second  partial  cycle  and  complete  cycle  of  motion  for  two-layer, 
cross-ply  plates  having  b/h=10  by  closed-form  and  finite-element 
methods* 


a/b 

aub^P/E^h3)5* 

^*2b2(P/Ef2h3)ls 

cob  2  (P/E 

f2h3)‘i 

C.F. 

F.E. 

C.F. 

F.E. 

C.F. 

F.E. 

Aramid-Rubber: 

0.5 

19.38 

20.23 

13.88 

14.55 

16.18 

16.93 

0.6 

14.65 

15.32 

11.05 

11.69 

12.60 

13.26 

0.7 

11.60 

12.17 

9.353 

9.807 

10.35 

10.86 

0.8 

9.537 

9.825 

8.269 

8.635 

8.860 

9.192 

0.9 

8.088 

8.488 

7.543 

7.879 

7.806 

8.172 

1.0 

7.038 

7.386 

7.038 

7.364 

7.038 

7.375 

1.2 

5.66T 

5.928 

6.402 

6.727 

6.008 

6.302 

1.4 

4.838 

5.045 

6.037 

6.356 

5.371 

5.625 

1.6 

4.313 

4.536 

5.812 

6.088 

4.951 

5.199 

1.8 

3.960 

4.116 

5.655 

5.910 

4.658 

4.852 

2.0 

3.712 

3-909 

5.551 

5.821 

4.449 

4.677 

Polyester-Rubber: 

0.5 

19.12 

19.81 

15.95 

16.61 

17.39 

18.07 

0.6 

14.42 

14.98 

12.26 

12.79 

13.25 

13.80 

0.7 

11.43 

11.92 

10.04 

10.45 

10.69 

11.14 

0.8 

9.435 

9.855 

8.632 

9.014 

9.016 

9.416 

0.9 

8.059 

8.421 

7.711 

8.144 

7.881 

8.280 

1.0 

7.084 

7.406 

7.085 

7.394 

7.085 

7.400 

1.2 

5.856 

6.111 

6.337 

6.613 

6.081 

6.352 

1.4 

5.164 

5.407 

5.928 

6.193 

5.520 

5.773 

1.6 

4.748 

4.986 

5.694 

5.928 

5.178 

5.416 

1.8 

4.485 

4.693 

'  5.543 

5.778 

4.958 

5.179 

2.0 

4.310 

4.518 

5.435 

5.688 

4.807 

5.036 

+Here  u>i  and  <o2  denote  the  frequencies  corresponding  to  the  first  and 
second  portions  of  a  cycle,  respectively,  and  u  denotes  the  effective  frequency 
for  an  entire  cycle. 


(a)  First  half  cycle 


(b)  Second  half  cycle 


Fig.  1  Bimodulus  action  during  the  two  half  cycles  of  motion  of  a  single¬ 
layer  bimodulus  plate.  Shaded  material  is  in  longitudinal  tension. 


(a)  First  portion  of  cycle 


(b)  Second  portion  of  cycl 


Fig.  2  Bimodulus  action  during  the  two  portions  of  motion  of  a  two-layer 
plate  In  the  fundamental  mode  of  vibration.  Bottom  layer  is  in 
x  direction  (0°),  top  layer  is  in  y  (90°).  Shaded  portions  are  in 
tension  in  the  respective  fiber  directions. 
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